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ABSTRACT 

A new harassment model is presented that models the complex, and dynamical tidal 
field of a Virgo like galaxy cluster. The model is applied to small, late-type dwarf 
disc galaxies (of substantially lower mass than in previous harassment simulations) as 
they infall into the cluster from the outskirts. These dwarf galaxies are only mildly 
affected by high speed tidal encounters with little or no observable consequences; 
typical stellar losses are < 10%, producing very low surface brightness streams {fj,B > 
31 mag arcsec"^), and a factor of two drop in dynamical mass-to-light ratio. Final 
stellar discs remain disc-like, and dominated by rotation although often with tidally 
induced spiral structure. By means of Monte-Carlo simulations, the statistically likely 
influences of harassment on infalling dwarf galaxies are determined. The effects of 
harassment are found to be highly dependent on the orbit of the galaxy within the 
cluster, such that newly accreted dwarf galaxies typically suffer only mild harassment. 
Strong tidal encounters, that can morphologically transform discs into spheroidals, 
are rare occurring in < 15% of dwarf galaxy infalls for typical orbits of sub-structure 
within ACDM cluster mass halos. For orbits with small apocentric distances (<250 
kpc), harassment is significantly stronger resulting in complete disruption or heavy 
mass loss (> 90% dark matter and > 50% stellar), however, such orbits are expected 
to be highly improbable for newly infalling galaxies due to the deep potential well of 
the cluster. 

Key words: methods: numerical — galaxies: dwarf — galaxies: evolution — galaxies: 
haloes — galaxies: interactions — galaxies: kinematics and dynamics 



1 INTRODUCTION 

Dwarf galaxies are easily the most common type of galaxy 
in the Universe, and yet their very origin remains a sub- 
ject of much debate. They can be found in vast numbers 
in the cluster environment - our nearest large cluster, the 
Virgo cluster, contains over a thousand dwarfs. The under- 
standing of the formation and evolution of cluster dwarfs is 
hampered by their intrinsic properties of low luminosity and 
low surface brightness causing them to be challenging ob- 
jects to observe, and study. Intuitively, it might be expected 
that their low mass, and hence shallow potential wells, cause 
them to be extremely sensitive to their local environment. 
As such, cluster dwarfs aught to provide sensitive indicators 
to the effects of mechanisms occurring in the high density 

environment. 

The Virgo cluster catalogue (jBinggeli et al.lll985l \ con- 
tinues to provide an invaluable tool for the study of dwarf 
galaxies in the Virgo cluster. Complete down to appar- 
ent magnitude Bt >18, and with a suprisingly accurate 
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morphological classification that has, on the whole, stood 
the test of time, the Virgo cluster catalogue (VCC) clearly 
demonstrates that dwarf galaxies are far from a homogenous 
family of galaxies. Dwarf ellipticals are the most common 
type of cluster dwarf - typically containing no current star 
formation, and an old stellar population. Dwarf ellipticals 
can be additionally split into further categories, with some 
containing a compact, central stellar nucleus - the nucle- 
ated dwarf ellipticals, and dwarf spheroidals - generally con- 
sidered to be low luminosity analogues of dwarf ellipticals. 
Not all the dwarfs are 'red and dead' with many containing 
active star formation, with knotty HII regions clearly vis- 
ible in the Las Campanas observatory photographic plates 
utilised to create the VCC. These star forming dwarfs are 
further split into sub-categories; dwarf irregulars containing 
atom ic gas, and with evidence of rotation (jvan Zee et aU 
and blue compact dwarfs whose gas and star forma- 
tion i s centrally concentra ted and also show significant rota- 
tion (|van Zee et al.|[200lh . A further classification is made 
for 'transition- types', objects whose properties are some- 
where mid-way between dwarf ellipticals and dwarf irreg- 
ulars, and may form a possible evolutionary link between 
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the two classes. The large variety in dwarf galaxy properties 
makes a single theory for their creation highly challenging. 

With the release of SDSS data set, covering the Virgo 
region, further progress was made in the observation of clus- 
ter dwarfs. A series of pap ers (iLisker et al.ll2006llLisker et al.l 
l|2006l ) lLisker et al.l (|2007l )) have provided a quantitive study 
of over 400 Virgo dwarf galaxies, suggesting that the VCC 
dwarf classifications can be further sub-categorised. The im- 
proved imaging quality, and additional colour information 
have shown that many dwarf ellipticals thought to be typ- 
ically 'red and dead', infact contain blue centres indicating 
recent or current central star-formation. Furthermore, dwarf 
ellipticals that appear smooth and featureless, present intri- 
cate spiral patte rns and disks when an un-sharp mask tech- 
nique is applied. iLisker et al.l (|2007h stress the importance of 
recognising these sub-categories and considering them sep- 
arately when attempting to un-ravel dwarf galaxy origins, 
and conclude that multiple environmental mechanisms are 
required to explain the diverse properties of cluster dwarfs. 

Current suggested dwarf galaxy formation and evo- 
lution scenarios are numerous and varied. The ACDM 
paradigm suggests that current day dwarf galaxies are 
the left over debris from hierarchical merging of haloes. 
While ACDM simulations predicts the presence of an ex- 
tended dark matter halo about each dwarf galaxy, they 
make less detailed predictions about properties of the bary- 
onic content of a galaxy. It is necessary to include phys- 
ical mechanisms that influence the visible components of 
dwarf galaxies in order to bring the simulations into agree- 
ment with observations (Bullock et al. 2000). These phys- 
ical mechanisms can be broadly split into two categories; 
global mechanisms that affect all dwarf galax ies regardless 
of th eir location, such as supernovae feedback (jPekel fc Silkl 
Il986l ). and environmentally dependent mechanisms. One 
such environmental mechanism, known as harassment, is 
believed to be of importance in high density environ ments 
ijMoore. Lake fc KatdflOoi . iMastropietro et al]|2005h . Nu- 
merous high speed tidal encounters within a cluster act to 
morphologically transform medium mass, low surface bright- 
ness disc galaxies into smaller galaxies - cluster dwarfs. A 
second scenario, ram pressure stripping, acts on in-falling 
galaxies that are already dwarfs - late-type dwarfs - and 
that this c an produce ga l axies that resemble cluster dwarf 
ellipticals l|van Zee et all |2004 iBoselh et all 120081 ). These 
two scenarios have one thing in common - a newly in- 
falling disc galaxy is subjected to cluster environmental 
mechanisms causing transformation into an object resem- 
bling a cluster dwarf elliptical. Indeed, radial velocity mea- 
surements of cluster dwarf ellipticals in the Virgo cluster 
do indicate that they are far from a virialised population 
within the cluster potential, be aring signatures of a re- 
cent in-fall (jConselice et al.ll200ll ). Meanwhile dwarf galax- 
ies obey their own morphology-density relationship; early 
-type dwarfs have been said to be the most strongly clus- 
tered of all galaxi es, whereas gas-rich d warfs are the most 
weakly clustered l|Ferguson fc Binggelil (|l994l )). As cluster- 
centric distance increases, the ratio of dwarf e llipticals to 
dwarf irregulars decreases (jSabatini et al.ll2005h . These two 
observations combined may suggest that many cluster dwarf 
ellipticals have their origin in infalling late-type dwarfs. Fur- 
ther evidence for galaxy transformation of discs is present in 
the form of hidden disc features, in some cases even includ- 



ing fi ne spiral arm structure l|Lisker et al.ll200"^ . ILisker et al.l 
[2003) in Virgo cluster dwarf ellipticals. It is as if the orig- 
inal discy nature of the galaxies has yet to be completely 
removed by the cluster environmental mechanisms. This is 
additionally supported by evidence of significant rotation i n 
a sub-sample of the dwarf ellipticals in Ivan Zee et al.l(|2004 ). 
Remarkably the rotating dwarf ellipticals, and dwarf irregu- 
lars, both appear to obey the TuUy-Fisher relation for bright 
late-type spirals. If we assume that a substantial fraction of 
present-day dwarf ellipticals are indeed accreted as dwarf 
discs into the cluster, the next question is whether it is ha- 
rassment or ram pressure stripping that is the key environ- 
mental influence driving the tranformation. As ram pressure 
stripping has been shown to remove significant amounts of 
atomic hydrogen from giant spiral galaxies wit hin the Virgo 
cluster (|Chung et al.ll2008l . rVoilmer et al.ll2002l '). it should be 
of little surprise that ram pressure can strip small late-type 
dwarfs of their HI gas content, causing a cessation of star- 
formation, and evolving the dwarf s from blue to red c olours, 
typical of cluster dwarf ellipticals l|Boselli et aLlfiooii ). How- 
ever, some form of physical morphological transformation 
may additionally be required to transform dwarfs from thin 
discs to spheroidals. One suggestion is that both mechanisms 
are required - ram pressure halts the star formation, while 
harassment thickens up the disc. 

In the literature l|Moore et al.lll99^ . harassment is the 
effects of repeated and numerous long range tidal encoun- 
ters in the potential well of the cluster. The tidal forces ex- 
perienced by a galaxy within a cluster can be broadly split 
into two categories - tidal heating, and tidal shocking. The 
key difference between these categories is the time-scale on 
which the galaxy experiences the tidal force. For tidal shock- 
ing the time-scale is short - for example in a high speed tidal 
encounter between two cluster galaxies, a strong but short- 
lived tidal force may be experienced by the two galaxies. 
Tidal heating occurs when the time-scale is far longer - for 
example a cluster galaxy that passes close to the cluster cen- 
tre spends a significant length of time under the influence of 
the deep potential well of the cluster, and may be entirely 
dismantled if the resulting tidal forces are in excess of the 
galaxy's self-gravity. Harassment is the combined effects of 
both tidal heating (from the tidal forces associated with the 
deep potential well of the cluster) and tidal shocking (from 
the tidal forces associated with the high speed galaxy-galaxy 
encounters) on a cluster galaxy. 

The aim of this work is to study the effects of harass- 
ment on newly infalling, late-type dwarf irregulars. This is 
accomplished using a new, and fast algorithm to model a 
galaxy clusters tidal fields that produces similar results to far 
more complex models set in a cosmological context. In con- 
cordance with cosmological models, both dwarf gala:xy mod- 
els and the cluster tidal field utilise cuspy cold dark matter 
haloes. The dwarf galaxy models are significantly lower mass 
(~10^ - 1 X 10^° Mq) than previous studies of harassment 
on dwarf galaxies. IMastropietro et"aL I l|2005l ) dwarf models 
were of mass 7 x 10^° Mq in comparison. Studying the effect 
of harassment on dwarfs in a lower mass regime is interesting 
for a number of reasons - obviously it extends harassment 
studies to cover a wider range of the parameter space. Addi- 
tionally while larger dwarfs may b e harassed into smaller ob - 
jects resembling dwarf ellipticals (|Mastropietro et al.ll2005l ). 
the fate of smaller mass dwarfs is less clear. A weakly ha- 
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rassed low mass dwarf may be able to complete a morpholog- 
ical transformation into an object resembling a dwarf ellip- 
tical without significant mass-loss. Alternatively, a heavily 
harassed low mass dwarf may be entirely destroyed in the 
process. 

For these low mass dwarfs, the significance of tidal heat- 
ing from the potential well of the cluster alone can be esti- 
mated using a Roche limit analysis. When the tidal forces 
of the cluster potential well exceeds the dwarf galaxy's own 
self-gravity, the dwarf will be entirely dismantled. For our 
standard dwarf model, and for the cluster potential utilised 
(see Section 13)), this can be expected to occur at a radius 
of less than ~ 50 kpc within the cluster. Hence only dwarfs 
on extremely plunging orbits are significantly infiuenced by 
tidal heating. Hence for the orbits we consider (see Section 
[5| whose peri-cluster distance is ~ 200 kpc, the effects of 
tidal heating alone are expected to be weak. 

We can also analytically estimate the effects of tidal 
shocking from galaxy-galaxy encounters for our low-mass 
dwarf models. Although it seems intuitive, that the lowest 
mass galaxies might be expected to be the most significantly 
affected by tidal encounters, this is not necessarily the case. 
The strength of the effects of a high-speed encounter can 
be expected to depend on the dynamical time at the scale- 
radius, Ts of the dwarf {tdyn ~ \{ri/GM{< rs)\^ , where 
G is the gravitational constant. When this dynamical time 
is long compared to the passage-time of a passing harasser 
gala:xy, then the the dwarf can respond effectively to the 
encounter. It is interesting to note that in the dwarf mod- 
els considered here, for material surrounding the disk of the 
galaxy, tdyn is of the same order as the passing time of a 
harasser galaxy. Despite the lowered mass, they are signif- 
icantly more compact, and hence resilient to harassment. 
These properties are not peculiar to our models - in this 
mass regime, this is to be expected. 

Although tidal stripping and shocking are expected to 
be weak for the dwarfs we consider, strong tidal shocking 
between galaxies in combination with tidal heating can still 
produce significant effects. These effects are studied in de- 
tail on model dwarf galaxies in Section [S] We attempt to 
quantify the likelihood of significant harassment in section 
(6] and test how the likelihood varies for a variety of orbits 
in Section [6.41 

The paper is organised as follows. In section 2, the numerical 
code utilised in all simulations is discussed, and in Section 3 
the methods used to build model dwarf galaxies are shown. 
The harassment model is presented in section 4. In section 

5, the infiuence of harassment on 11 in-falling dwarf galaxy 
models are shown. The statistical infiuences of harassment 
are found by means of a Monte-Carlo simulation in section 

6, and tested for a variety of orbits. 



2 THE CODE 

In this study we make use of 'gf l|Williams fc NelsonI 
l20QlllWilliamj 1199 8'). which is a Treecode-SPH algorithm 
that operates primarily using the techniques described in 
iHernauist fc Kat j (|l989l '). ' gf ' has been parallelised to oper- 
ate simultaneously on multiple processors to decrease simu- 
lation run-times. While the Treecode allows for rapid calcu- 
lation of gravitational accelerations, the SPH code allows us 



to include a HI gas component to our dwarf galaxy models. 
In all simulations, the gravitational softening length, e, is 
fixed for all particles at a value of 100 pc, in comm o n wit h 
the harassment simulations of iMastropietro et al.l (|2005l ) . 
Gravitational accelerations are evaluated to quadropole or- 
der, using an opening angle 9^ = 0.7. A second order in- 
dividual particle timestep scheme was utilised to improve 
efficiency following the methodology of iHernguist fc Kata 
(1989). Each particle was assigned a time-step that is a 
power of two division of the simulation block timestep, with 
a minimum timestep of ^0.5 yrs. Assignment of time -steps 
for co Uisionless particles is controlled by the criteria of lKat j 
|l993), whereas SPH particle timesteps are assigned us- 
ing the minimum of the gravitational time-step and the 
SPH Courant conditions with a Courant c onstant, C=0.1 
(Her riauist fc Kat3ll989l ). As discussed in IWilliams et aTl 
(200^, the kernel radius h of each SPH paricle was allowed 
to vary such that at all times it maintains between 30 and 
40 neighbours within 2h. In order to realistically simulate 
shocks with in the SPH model, th e artifi cial viscosity pre- 
scription of iGingold fc MonaghanI (|l983t ) is used with vis- 
cosity parameters (a,/3) = (1,2). The equation of state for 
the gas component of the galaxies is isothermal with a bulk 
velocity dispersion of 7.5 km s~^, based on the measured ve- 
locity dis persion of molecular c louds in the local interstellar 
medium (I St ark fc Brand [l98^ ). By choosing an isothermal 
equation of state, we are intrisically assuming that stellar 
feedback processes are balanced by radiative cooling produc- 
ing a constant velocity dispersion, 'gf contains a simple star 
formation prescription, where each SPH particle converts it's 
mass into stellar mass at a rate controlled by a Schmidt law 
with a power of 1.5. When the global stellar mass formed 
in this manner reaches a specified mass, a specified number 
of new N-body star particles are formed. The new particles 
have their position chosen statistically such that recent star 
formation history of individual SPH particles is t aken into 
accou nt. For details of code testing, please refer to I William j 

(jliii). 



3 BUILDING MODEL DWARF GALAXIES 

Our dwarf galaxie s models consist of 3 component s; an NFW 
dark matter halo (|Navarro. Frenk fc WhitellTgg^ ). an expo- 
nential disc of gas and one of stars. The methods used to 
form each component will be discussed in the following sec- 
tions. 

3.1 The dark matter halo 

In the current standard model of galaxy formation galaxy 
mass is dominated by an unseen dark matter component, 
surrounding the disc. This provides the additional grav- 
itational force required to explain the observed velocities 
while s imultaneously produc ing their observed flat rotation 
curves. iNavarro et "alT(ll996l ') suggests that the density pro- 
files of dark haloes have a universal shape (the NFW profile) , 
for haloes that range in mass from dwarf galaxies to galaxy 
clusters. 

The NFW profile has the form: 

pNFw{r) = j_Y (1) 
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where Va is a characteristic radial scale-length. The profile is 
truncated at the virial radius, r2oo = VsC. Her e c is the con- 
centration parameter (jLokas fc MamonlboOll ;). c is found to 
have a range of values in cosmological simulation, however 
there is a general trend for higher val ues in less mass i ve sys - 
tems with some scatter - see figure S. lNavarro et al.l l|l996l ). 
Typical values for cluster-mass objects is c ~4, whereas 
dwarf galaxies can have c ~ 20. For our standard dwarf 
galaxy model we choose a total mass of 1O^"M0, and c=20. 
Positions and velocities are assigned the dark ma tter par- 
ticles using the distribution function described in IWidrowl 
(|2000l ). 

Dark matter haloes produced in this manner are tested 
by evolving the initial conditions using the numerical sim- 
ulation code, and the evolution of the density profile is ob- 
served. Tests show that after 0.5 Gyrs, transient effects have 
settled, and the halo remains stable. Due to the abrupt trun- 
cation of the density distribution at r2oo, ~ 2% of halo mass 
can expand outwards at the truncation, increasing it's radii 
by ~ 25% but this causes negligible effect on the density 
distribution surrounding the disc. 



3.2 The galaxy discs 

Within the centre of the NFW halo is a rotating stellar and 
gas disc. The stellar and gas disc of typical spiral galaxies 
are observed to be exponential, albeit with roughly a fac- 
tor of two difference betw een the gas and stellar scalelength 
l|Boselh fc Gavazzill2006l ): 



E(r) = T.oexp{r/ra 



(2) 



where E is the surface density, Eo is central surface density, 
r is radius within the disc, and is the scale-length of the 
disc. 

In order to choose a reasonable scale-length fo r the stel- 
lar co mponent of the disc, we follow the recipe oflMo et al.l 
l|l998 l) for producing discs within the ACDM paradigm. This 
recipe has been shown to reproduce the TuUy-Fisher rela- 
tionship's slope and scatter, as well as the observed scatter 
in the size-rotation velocity plane. Having first chosen a halo 
mass, the disc mass Md is assumed to be a fixed fraction 
of the halo mass. Additionally, the angular momentum of the 
disc is assu med to be a fixed f raction jd of the halo's angular 
momentum. IMo et al] (liggsl ) find nid ~ jd and md ^0.05 
to reproduce observations, and we choose mti=0.05 for all 
gala:xy models. The scale-length of the disc is then fully de- 
fined by the properties of it's host halo: 



2^200 \E\^/^ 



rud 



(3) 



where M200, V2oo,and E are total mass, circular velocity, 
and energy of the halo, and A is the spin parameter of 
the halo. We choose A=0.05 for all disc models presented 
in this paper u nless otherwise stated (the mean A in the 
I Mo et al] (|l998l ) probability function) . For our standard 
galaxy model, this produces an rs=0.86 Kpc. Knowing this, 
we can solve for E(r) in Equation [5] using Eo = Md/^Tirj,. 
The disc scalelength for the gas disc is then simply the stel- 
lar disc scalength multiplied by a factor of two. In order to 



set up the positions of stellar and gas particles, we integrate 
over E(r-) in concentric rings using a standard integral such 
that 



M(r') 



T,(r)2TTrdr = 2nT,ors 



-r 



+ 1 

(4)' 

As with the dark matter halo, we can solve for r'(A/) nu- 
merically, then by producing a random number between 
and 1 for the normalised M, we are provided with the ra- 
dius of each particle. These particles are initially laid down 
in a single plane. Real discs have a finite thickness, but we 
cannot solve for this until we have defined the velocity dis- 
persion throughout the disc. This is n ecessary in or der to 
ensure that the disc is Toomre stable l|Toomre|[l96i ). The 
Toomre stability criterion is defined as 



Q 



3.36GE 



> 1 



(5) 



Here E is the surface density, ar is the radial velocity dis- 
persion, and K is the epicyclic frequency defined, using the 
epicyclic approximation (Springcl fcJWhitc 1999). Next we 

use = ^ where 7^ = ^ (also ISpringel fc White! 
'1999I), and 0^ ^0.6-(j>r (|Shlosina!ii fc Noguchil (|l993t )). This 
completely defines the minimum values of the necessary ve- 
locity dispersions throughout the disc. In practice, Q > 1.5 
is required throughout the stellar disc to ensure stability. 
The gas disc has an intrinsic velocity dispersion, due to it's 
isothermal nature, that automatically sa tisfies the Toomre 
criteria at all radii. Once more following ISpringel fc WhiQ 

(|l999t ). we now use Zd = —7% for the vertical scale height 
of the disc, and distribute the particles vertically out of the 
disc following the Spitzer's isothermal sheet solution 



p{R,z) = ^^sech^(2/2d) 



(6) 



Finally, the circular velocities of disc particles are calcu- 
lated. The value of the potential is calculated in thin spher- 
ical shells for the combined dark matter, stars and gas in 
a one-off N-body calculation, including the effects of grav- 
itational softening. This is necessary as the discs own self- 
gravity can influence the radial orbits of disc particles close 
to the centre of the halo. The gradient of the potential is 
then used to calculate the circular velocity at each radius. 
As a result inner disc particles have raised circular velocities 
beyond that of the circular velocity of the halo alone at that 
radius. 

In all standard dwarf galaxy model simulations, a gas- 
rich dwarf irregular is simulated, containing 4 x 10* Mq 
of gas and 1 x 10* Mq of stars. Such gas rich systems 
can be found in Virgo cluster dwarf galaxy observations 
l|BoseUi fc Gavazzilbood ). However the paper's conclusions 
are additionally tested using a pure stellar galaxy disc model 
to verify the significance of the gas-to-stellar ratio on the in- 
fiuences of harassment - see Section 15.51 

To ensure long term stability, the complete three- 
component galaxy models are evolved in isolation for 2.5 
Gyrs. Transient effects are found to fully settle in <1 Gyr 
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Component 


Parameter 


Parameter value 


Halo 








Concentration c 


20 




Virial Radius r2oo 


44 kpc 




Particle number 


200000 


Total disc 


Mass 


5 XlO* Mq 




Particle number 


50000 


Stellar disc 


disc fraction 


0.2 




Particle number 


10000 




Scalelength r^j 


0.87 kpc 


Gas disc 


disc fraction 


0.8 




Particle number 


40000 




Scalelength 


1.74 kpc 



Table 1. Summary of key parameters of standard dwarf galaxy 
model 



and only evolved models are introduced into the cluster ha- 
rassment models. 

A summary of the standard dwarf galaxy models prop- 
erties is provided in Table [T] 



4 THE HARASSMENT MODEL 

The potential well from the gravitational field of a large 
galaxy cluster, such as Virgo, applies a significant force on 
nearby galaxies. The local group itself is thought to be in- 
falling towards the c luster at ~200 km s~^ as a r esult of the 
Virgo over-density l|Tammann fc Sandagd [l985l ). However, 
it is the tidal force of the cluster that will strip or deform an 
in-falling galaxy. Observationally, clusters have typical mass- 
to-light ratios of several hundred dSheldon et al.. ,2009 ). As 
a result, the combined mass of individual galaxies, even in- 
cluding their dark matter haloes, is not sufficient to total 
the mass of a cluster. 

As a result, we model the tidal field of a Virgo- like clus- 
ter using a two component model; the first is a dynamical 
tidal field associated with individual 'harasser' galaxies, and 
the second is a static tidal field to represent the rest of the 
mass of the cluster that is not associated with galaxies - 
referred to as the background cluster potential herein. The 
potential field of both the harassers and the background is 
represented using the analytical form for an NFW density 
distribution (|Lokas fc MamonI (|200ll 'l') 



ln(l + {r/r,)) 



(7) 



where Qc = l/[/7i(l + c) ~ c/(l + c)]. 

Modelling the total cluster potential as two compo- 
nents has an additional advantage - it is easy to separate 
and study the effects of the background cluster potential 
from the additional effects of high-speed galaxy-galaxy en- 
counters. The background cluster potential is fully defined 
once its virial mass and concentration have been chosen. We 
choose c=4, and M2oo = 1.6x IO^^Mq. This choice of cluster 
m ass is in agreement wi th the model Virgo cluster potential 
of IVollmer etal] l|200ll ) and with ga l axy li ne-of-sight ve- 
locity measurements in IBinggeli et al.l (119851 ) (~ 10^*Mq). 
Now each harasser galaxy has it's own individual NFW 
potential field. To make this field dynamical, a one-off N- 



body simulation of a cluster-mass halo is conducted, and 
the time-dependent coordinates of a fraction of the particles 
are logged. This log now provides the position-evolution of 
the centre of each of the harassing haloes, where each in- 
dividual galaxy potential field will be super-imposed. Once 
more we require a mass and concentration for each harasser's 
potential fiel d to be defined. For mass, we use the Schec- 
tor function llSchechteil 11976* with parameters provided in 
ISandage et al.1 (|l985l ). and assume a constant mass-to-light 
ratio, i.e. 



Mtot 



exp 



L_ 



dL (8) 



where a = —1.25, and (M/L)k = 20. The fitting param- 
eters provided in |Sandago ot al.' ("1985) are based on the 
Virgo cluster catalogue (Binggeli ct al. 1985;), and are there- 
fore complete down to a minimum galaxy luminosity of ~ 
2.5 X 10^ Lq (assuming a distance modulus of 31.0 for Virgo). 
We choose to resolve all harasser galaxies down to one-tenth 
the mass of our infalling dwarf model (IO^Mq). Completing 
the integral in Equation |8] produces a total mass in galaxies 
that is 14.5% of the total cluster mass. This is in reasonable 
agreement with the ~ 10% percentage of a cluster's mass 
loc ked up in s ub-str ucture, found in ACDM simulations - 
see iGill et al.l \2004 ]. This process produces 733 harassing 
galaxies in total. The cummulative distribution of galaxy 
masses produced can be seen in Figure [l] The vast major- 
ity of harasser galaxies are dwarfs, with ~ 65% of the total 
number of harassers having a mass less than the standard 
dwarf galaxy model, while only ~ 9% of the harassers have 
a mass greater than 10 times that of the model dwarf. The 
maximum mass galaxies produced are ~ lO^^Af0. The value 
of the mass-to-light ratio chosen is consisten t, if not an up- 
per l imit for galaxies in this mass-range (see lGilmore et ahl 
l2007h . The concentration c of each halo is chosen to follow 
the trend of higher concentration for lower mass objects as 
found in ACDM simulations. This is a chieved by a fit t o con - 
centration values with mass found in iNavaxro et al.1 (|1996D 
producing c — —S\og{Mhaio{MQ)) + 52. This fit produces 
c ~ 15 for a 6 X lO" Mq mass halo, and c ~ 20 for a 1 x 10^" 
Mq mass halo. 

It should be noted that the use of analytical po- 
tentials for the harassing galaxies tidal fields has advan- 
tages and disadvan t ages o ver p revious haras s ment models. 
iMastropietro et al.1 (j^ol) and iMoore et all l|l999l ) utihse 
full ACDM sim ulations of t h e dar k matter in cluster-mass 
objects, while , Moore et al.l (|l998l ) utilises softened point 
masses to represent the tidal fields of harassers. These previ- 
ous methods have the advantage that the harassing galaxies 
can interact and react to the tidal potential of the harassed 
galaxy model, whereas analytical potentials on fixed paths 
of motion do not. This essentially violates Newton's third 
law - the dwarf can respond to the cluster potential field, 
but the cluster potential field cannot respond to the dwarf. 
However, primarily we are interested in studying the effects 
of the harassers on the model galaxy and not vice versa. As 
will be shown in Section [6.31 encounters velocities are typ- 
ically very high, and as a result the tidal forces in a high 
speed encounter are brief, and strongly peaked. The abil- 
ity for the harasser to alter it's trajectory in response to 
the dwarf models tidal field is inconsequential to the brief. 
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Figure 2. A tidal history for Run 4. Tidal forces are measured 
in units of the isolated dwarf galaxy model's Roche limit at the 
edge of it's stellar disc. The Gaussian-like curve represents the 
contribution from the back ground cluster potential, and points 
are the contribution from individual galaxy-galaxy encounters. 



strong tidal forces that the model is subjected to at these 
velocities. Typically the dwarf model suffers a slight devia- 
tion in it's orbit as the result of a strong tidal encounter, and 
the same might be expected for the harasser galaxy. How- 
ever this will have a negligible effect on the strength, and 
shape of the peaked tidal forces due to the rapid passing 
time, and high velocity of the encounter. Additionally, the 
use of analytical potentials has advantages - the tidal forces 
have effectively infinitely small spatial resolution. Only the 
self-gravity of the galaxy model's own particles have a lim- 
ited spatial resolution due to softened gravity. Analytical po- 
tentials are highly computationally cheap, allowing a large 
number of simulations to be conducted rapidly as is a neces- 
sity for harassment simulations due to their highly stochastic 
behaviour (see Section |6}. 



5 EFFECTS OF HARASSMENT ON MODEL 
DWARF GALAXIES 

The harassment model is utilised for 8 infalls (Runs 1-8) of 
the standard dwarf galaxy model. Initially, the dwarf galaxy 
models are positioned at the outskirts of the cluster, and 
randomly orientated on a shell at the virial radius (rapo=l-l 
Mpc), such that each dwarf galaxy is presented with a differ- 
ent set of tidal encounters as it infalls. They are each given 
an azimuthal velocity of 250 km s~^. If all harassing galaxies 
are removed, this produces a plunging orbit with a pericen- 
ter of ~ 200 kpc and a peak velocity of ~ 1600 km s~^. This 
choice of orbit is consistent with the rang e of orbits utilised 
for infalling galaxies in ,Vollrner e t al. (200^) . The orbit is 
chosen to model a dwarf galaxy infalling into the cluster for 
the first time at the current epoch. Models are evolved for 
2.5 Gyrs providing time for the infalling galaxy to make one- 
pass of the cluster centre on this orbit. In order to complete 
a second pass, the model would need to be evolved for > 6 
Gyrs in which case the assumption of constant cluster mass 
becomes invalid. Hence the model is most applicable to a 
first infall into the cluster environment. The effects of less 
plunging orbits are simulated and discussed in Section [6] 



5.1 Tidal histories 

As we know both the position-evolution and the proper- 
ties of all harasser galaxies, we can measure the exact tidal 
force, from the background cluster potential and from ha- 
rasser galaxies, felt by an infalling galaxy at any instant. 
All tidal force measurements are presented in units of the 
Roche limit of the isolated dwarf galaxy model at the edge 
of its stellar disc - a discussion of the use of Roche limits 
is provided in the following section. In Figure (2] the tidal 
evolution of Run 4 is presented as a function of time. We 
refer to these figures as 'tidal histories', and Run 4 can be 
considered a typical infall for reasons that will be further 
discussed in Section [S] Clearly the background cluster po- 
tential alone doesn't produce tidal forces beyond the dwarf 
galaxy's Roche limit, even at the peri-center of the orbit. 
As discussed and quantified in the Introduction, this would 
require a far more deeply plunging orbit. However, it is the 
sum of both the background and the individual encounters 
that can produce tidal forces exceeding the Roche limit. 

5.2 Mass loss 

Although the net tidal force rises above the Roche limit, it 
does so only briefly and hence the complete dismantling of 
the dwarf does not occur. This form of tidal force is therefore 
in the 'tidal-shocking category' (see Introduction for further 
discussion). The use of tidal forces expressed in units of the 
Roche limit of the dwarf is directly relevant for tidal heat- 
ing, but not so for tidal shocking. As a result, it should be 
clearly noted that the presenting of all tidal forces in units of 
the Roche limit serves merely as a more intuitive and com- 
prehendable scale or yard-stick for the strength of the high 
speed tidal encounters. In figure [S] the number evolution of 
dark matter, and stellar particles within a sphere of radius 15 
kpc centered on the stellar disc is presented. Clearly the stel- 
lar disc is far from completely dismantled, and loses < 10% 
of it's original mass. The dark matter halo is affected far 
more significantly, losing ~ 60% of it's original mass. To test 
that this quantification of dark matter loss is robust to in- 
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Figure 3. Number evolution of dark matter particles (solid), and 
star particles (dashed), within a sphere of radius 15 kpc centred 
on the stellar disc for Run 4. The dotted line shows the final 
number of dark matter particles at the end of a simulation with 
the background cluster potential alone, i.e. no harassers included 

creased resolution, we repeat the infall witli tlie gravitational 
resolution improved by a factor of 10 (the softening length, 
e is reduced by a factor of 10 to 10 pc). We also improve the 
dark matter mass resolution by a factor of 2.5 by conducting 
the high resolution simulation with a dark matter halo con- 
sisting of 500,000 dark matter particles. The high resolution 
models are found to have negligible differences in dark mat- 
ter and stellar mass losses from the lower resolution studies. 
To confirm how much of the dark matter loss is from the 
background potential alone, and how much is lost because 
of the inclusion of the harassing galaxies, the simulation is 
repeated without the harassing galaxies. In this simulation, 
the galaxy models loses only ~ 25% of its dark matter - 
shown as dotted line in Figure O Therefore galaxy-galaxy 
interactions are responsible fo r more than half of the dark 
matter losses. As discussed in lMoore et al.l l|l998h . the rea- 
son that the stars are so less significantly affected is due to 
their circular, centralised orbits whereas the eccentric orbits 
of dark matter particles render them far more susceptible 
to tidal disruption. The net result of enhanced dark matter 
losses over stellar losses is a decrease in dynamical mass-to- 
light ratio by roughly a factor of 2. With currently available 
data sets, it would be a significant observational challenge 
to detect a systematic difi^erence in dynamical mass-to-light 
ratio of isolated dwarfs in comparison to cluster dwarfs at 
this level, due to the inherent observational and theoretical 
difficulties in obtaining these ratios. Such a test would pro- 
vide a significa nt test of our under standing of dark matter. 
As discussed in lMoore et al.l (| 19981 ). it is difficult to see how 
cluster galaxies can avoid losing a sizeable fraction of their 
dark matter, within the ACDM paradigm. 

5.3 Little morphological transformation 

The final harassed stellar disc does not appear to be sig- 
nificantly transformed from the original disc. The discs re- 
mains flattened and dominated by rotation. The isolated 
dwarf galaxy has a peak circular velocity Vpeafc '^SO km 
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Figure 4. B-band surface brightness plot of face-on stellar disc 
in Run 3 after 2.5 Gyrs of harassment. The grey-scale colour bar 
indicates surface brightness in units of mag arcsec"^. The box 
size is 5 X 5 kpc. White isochrones are displayed with values of 
25, 28, and 31 mag arcsec"^ (from galaxy centre to edge of disc) 

s~^, and a velocity dispersion a out of the plane of the disc 
~ 10 km s~^. After harassment, Vpeak falls slightly, and a 
is raised to 14 km s~^. This has a net effect of reducing the 
(V^eaft/o") ratio by less than half to 4.6, therefore remain- 
ing significantly dominated by rotation. 

The post-harassment stellar discs, which show induced 
spiral structure and bars, were initially smooth and fea- 
tureless - see Fig ure [J] It is interesting to note that in 
iLisker et al.l l|2007l ) , spiral structure in dwarf ellipticals is ob- 
served, and presented as evidence for the origin of dwarf el- 
lipticals in disc galaxies. However, in these simulations such 
features originate because of high speed tidal encounters. The 
stars that are stripped are arranged along very low surface 
brightness streams, that move along similar orbits to that 
of the galaxy though the cluster. Note that in Figure |4j the 
extended low surface brightness streams are not visible even 
at surface brightness limits of /is > 31 mag arcsec"^. To 
form these image s, a stellar mass-to-light ratio of ~ 6 is 
assumed following iMastropietro et al.l (|2005l ). These surface 
brightnesses are well beneath the surface brightness limits of 
current optical Virgo galaxy surveys. For exa mple, the INT 
wide- field survey has hb ^ 2& mag arcsec"^ (iDavi es et al.l 
l2005h . 

Although we have concentrated on one typical harass- 
ment infall simulation, 6 out of 8 of the dwarf galaxy infalls 
are affected in a very similar manner to the results presented 
above. In Table[2j the effects of harassment in terms of mass 
loss and change in (Vpeak /o") is summarised for all 8 simula- 
tions. The exception is Run 5 and Run 6 for which, in both 
cases, the galaxy model has at least one encounter of ^ 4 
times the Roche limit or greater. To better understand the 
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Run 


Tides 


DM 


Star 




{M/L)dj,„ 


1 


2.4,2.5,1.4 


30 


82 


45 


37 


2 


2.5 


59 


100 


59 


59 


3 


1.1,1.8,2.3 


65 


94 


69 


70 


4 


1.1,1.2,1.6 


41 


93 


55 


44 


5 


1.8,1.7,5.6 


5 


22 


30 


22 


6 


3.7,1.9,1.8 


8 


34 


44 


23 


7 


1.6,1.1,1.2 


33 


91 


57 


36 


8 


2.3 


47 


98 


58 


48 



Table 2. Summary of 8 dwarf galaxy infalls; column 1 is run 
description, column 2 is summary of all tidal encounters with 
tidal strength greater than the Roche limit of the isolated dwarf 
galaxy, column 3 and column 4 are percentage of dark matter and 
stellar mass at the end of the simulation, found within a 15 kpc 
radius sphere, column 5 and column 6 is final (Vpeafe/o"): ^^nd final 
(M/L)^j^„ as a percentage of the orignal 



frequency of such encounters we utilise Monte-Carlo simu- 
lations (see Section |6]). 



5.4 Effects on gas and star formation 

In Figure [S] the effects of high speed tidal encounters on the 
gas component of the infalling dwarf can be seen. Prior to 
the infall (left panel) the gas disc is a smooth exponential 
disc. Tidal encounters can induce spiral structure, which is 
stronger than is seen in the stellar disc due to the dissipative 
nature of the gas ( central panel) . The deformation of the disc 
from a circular shape causes shocking and barring of the gas 
disc, resulting in enhanced star formation within the spiral 
arms. Loss of angular momentum along the spiral arms can 
additionally cause radial inflow of the gas into the central 
regions of the dwarf irregular. However the spiral structure 
is only temporarily induced and, once the driving mecha- 
nism has stopped, the disc becomes smooth once more {right 
panel), leaving only a central knot of gas. The short-lived en- 
hancement of gas densities within spiral arms and the cen- 
tral region results in a total galaxy star formation rate that 
is bursty in nature (Figure [6]). The star formation rates of 
the dwarf irregular models decline steadily as they convert 
their gas reservoir into stars. However tidal encounters can 
cause both brief lowering (when the gas disc is stretched 
and lowered in density), and raising (when the gas disc is 
compressed) of the star formation rates. It should be noted 
that the simplicity by which star formation and associated 
feedback processes have been modelled here, provides only 
qualitative predictions for star formation rates. For example, 
a prescription for super-nova feedback, and gas cooling is not 
specifically included. Instead, it is assumed that energy in- 
put from stellar feedback is balanced by radiative cooling to 
produced an isothermal, single phase gas representing the 
diffuse, atomic gas content of the galaxy. The artificial vis- 
cosity of the gas component does enable it to shock, and this 
is important for spiral-arm formation, and radial gas infiow 
towards the galaxy centre. We do not include a prescription 
for gas-ionisation from background U V photons that m ay 
be important in galaxies of this mass (iMaver et al.l l|2006l) ). 
Hence, the star-formation rates seen in these simulations 
can only be considered qualitive at best, based purely on 
the density of the gas in the disk. 



However, it is likely that the missing gas physics is of 
little consequence as infalling dwarf galaxies will be most 
significantly infiuenced by ram pressure stripping in the clus- 
ter environment. The effects of ram pressure s tripping will 
be in vestigated in detail in a subsequent paper l|Smith et al.l 
l201(t ). For now, we restrict our analysis to the use of a simple 
analy tical prediction for t he efficiency of ram pressure strip- 
ping (|Gunn fc Gottlll972l ). Despite its simplicity, the Gunn 
and Gott formalism has been proven to provide reason- 
able first order predictions for the gas truncation radius in 
agreement with numerous full, numerical simulations of ram 
pressure stripping (Abadi ct al. 1999, VoUmcr ct al. 200^, 
Mayer et al1l2006i . Uachvm et al.ll200l iRoediger fc BruggenI 
20071 ) . We use the intra-cl uster medium radial density pro- 
file o f IVoUmer et al] (|2002l ). for the varying infall velocity of 
our infalling dwarfs, to predict the radius within the cluster 
when the dwarf's gas disc will be completely stripped. This 
method predicts that the dwarf will be completely stripped 
at a cluster radius of ~ 350 kpc, before it completes it's first 
pass of the cluster centre. Thus, it is unlikely that the major- 
ity of infalling dwarfs will still contain appreciable amounts 
of gas before they experience harassment. The typical loca- 
tion within the cluster that an infalling dwarf experiences 
harassment is discussed further in Section [6.31 Additionally 
this excludes the possibility that a central nucleus could 
form from the radial inflow of gas seen in the harassment 
simulations (although it should be noted that nuclei tend 
to be observed mo re frequently in brighter dwarf galaxies 
(jLisker et al.1 12007| ) that may be able to maintain a trun- 
cated HI disk after ram pressure stripping, and hence could 
allow formation of nuclei in this manner). 



5.5 Dependency on disc properties 

Our harassment model is ideal for testing how a galaxys re- 
sponse to harassment depends on specific properties such as 
gas fraction, surface brightness and mass. As will be shown 
in Section [S] a tidal history such as that of Run 4 can be 
considered a statistically typical infall. Now different galaxy 
models can be made to infall along the same trajectory, and 
be subjected to an identical tidal history. Using this method 
we test how different dwarf galaxy's respond to harassment; 
(i) Model A contains no gas, (ii) Model B is a more low 
surface brightness disc, and (iii) Model C is 6 times less 
massive. 



5. 5. 1 Model parameters 

Model A is similar to the standard dwarf model, except it's 
disc mass is purely stellar. Each gas particle is replaced with 
an identical mass stellar particle. Model B is also similar 
except that we utilise a spin parameter A = 0.1. This is 
double the value of the standard dwarf galaxy model and 
forces the exponential disc scalelength of the stellar and gas 
disc to also double, to = 1.72 kpc (stars) and 3.48 kpc 
(gas). As a result, model B is 25% lower in surface brightness 
at all radii within the disc than in Model A. Model C has 
a dark matter halo of 1.6 x 10^ Mq. As it's dark matter 
halo is less massive, it's c oncentration is in creased slightly 
to c = 24. Following the IMo et al.1 l|l998l ) recipe for this 
halo, produces a stellar disc with an exponential scale length 
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Figure 5. Column density plot of facc-on gas disc in Run 4; {left) before harassment, {centre) post tidal encounter, {right) at the end 
of the simulation. The grey-scale colour bar indicates column density in units of Mq pc"'^. The box size is 10 X 10 kpc 
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Figure 6. Bursty star formation rates of the 8 infalling dwarf 
irregular models. Star formation rates are shown in Mq yr~^ 



rd = 0.43 kpc. Once more, the scalelength of the gas disc is 
double that of the stellar disc. 



5.5.2 Response to typical harassment 

In Table [3] the key effects of harassment are summarised for 
ease of comparison. Particles are counted within a 15 kpc 
radius, to calculate dark matter and stellar losses in Model 
A. A sphere of double this radius was used for Model B, to 
match the fact that it's stellar disc had doubled in diameter. 
A sphere of only 7.5 kpc was used for Model C to account 
for it's smaller stellar disc. 

Firstly, comparing the standard model to Model A, 
there appears to be very little change in the effects of ha- 
rassment in terms of mass loss, mass-to-light ratio, or stellar 
dynamics. A visual inspection of the stellar disc also reveals 
little difference, with mild enhancement of spiral structure 
in both models. The inclusion of a gas component has made 
little difference. This suggests thats the effects of harass- 
ment will not change significantly if ram-pressure stripping 
had removed the gas from the gala:xy prior to harassment 
(see Section [5.41 for further discussion). However, a deeper 



Run 


DM 


Star 


{Vpeak/f^) 


(M/L)dy„ 


Standard 


41 


93 


55 


44 


Model A (star only) 


42 


91 


61 


46 


Model B (LSB) 


41 


75 


34 


55 


Model C (low mass) 


37 


97 


73 


38 



Table 3. Summary of effects of typical harassment on varying 
dwarf galaxy models; column 1 is run description, column 2 and 
column 3 are percentage of dark matter and stellar mass at the 
end of the simulation, coiiimn 4 a-nd column 5 is final (V^eafe/f), 
and final (M/L),^^^ as a percentage of the orignal 



study of the combined effects of harassment a nd ram pres- 
sure stripping is deferred to a future paper l|Smith et al.l 
I2OIOI ). 

The same cannot be said for the low surface brightness 
disc of Model B. The halo has suffered similar losses as the 
standard model - this is to be expected as their haloes are 
identical and suffering the same tidal interactions. But the 
stellar disc suffered additional ~ 15% losses. The (V^eafc/o") 
ratio is also reduced by an additional 10% from the standard 
model. The simplest explanation is that a larger stellar disc 
represents a larger target for harassing galaxies. In addition, 
the self gravity of the disc is reduced by a factor of 25%. This 
is in agreement with the simulations of iMoore et al.l (Il999l ) 
where a low surface brightness giant spiral suffers consider- 
ably more harassment than a high surface brightness spiral. 
The surface brightness (or disc scale-length for exponential 
discs) appears to be an additional parameter controlling a 
galaxys sensitivity to harassment, as it controls the cross- 
sectional area of a galaxy to harassers. 

Finally, we compare the low mass dwarf (model C) to 
the standard model. Despite model C having considerably 
less mass, its dark matter halo has not suffered substan- 
tially greater mass-loss. Furthermore, it's stellar disc losses 
are a smaller mass fraction than in the standard model. The 
(Vpeak/i^) ratio suffers only a minor reduction in comparison 
to the other models. This further supports the conclusion 
that the cross-sectional area of a galaxy is a strong param- 
eter controlling the response of a galaxy to harassment. 
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Figure 7. Cummulative frequency plot of peak strength tidal 
encounters experienced by infalling dwarf galaxies on plunging 
orbits 



6 A MONTE-CARLO SIMULATION OF 
HARASSMENT 

Column 2 of Table [2] summarises the tidal histories of each 
infalling dwarf galaxy model. There is a large range in 
strength and frequency of tidal encounters. This brings into 
question just how commonly strong encounters occur, or in 
fact just how typical our 'typical' simulation really is. To 
attempt to answer this question, and to gain additional in- 
sight into the principles of harassment, we conduct a 1000 
galaxy infall Monte-Carlo simulation. 

Clearly, conducting 1000 full N-body/SPH simulations 
of infalling dwarf galaxy models would be extremely time- 
consuming. Instead 1000 N-body particles with mass equal 
to the dwarf galaxy (and with the same assumed Roche 
limit) are allowed to infall into the cluster model. Initially 
each particle is randomly positioned on a spherical surface of 
radius 1.1 Mpc. They are given an initial velocity of 250 km 
s~^ and the direction of the velocity vector on the spherical 
surface is randomised. Hence the initial orbital parameters 
(rapo and Vazi) of each particle's orbit is identical to that 
used in the full numerical simulations. However, now the 
range of tidal histories that an infalling dwarf galaxy en- 
counters is sampled to a far greater degree of completeness 
than is possible with only 8 galaxy model infalls. The tidal 
history of each particle is then recorded allowing for a sta- 
tistical analysis, and by comparison with the tidal histories 
of the full simulations, conclusions can be drawn. 



6.1 Frequency of strong encounters 

To understand, statistically, what the typical tidal history of 
an infalling dwarf galaxy appears like, we log the peak tidal 
force that each infalling Monte-Carlo particle is subjected to 
throughout its orbit. In Figure [71 a cummulative frequency 
plot of the peak strength encounter is shown. It can be seen 
that half of the infalling dwarfs never experience a peak tidal 
force greater than twice the Roche limit of the dwarf model. 
In fact peak tidal forces in the range of 1-2 Roche limits 
are the most common experienced by the infalling galaxies. 
In the full numerical simulations, both Run 4 and Run 7 
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Figure 8. Peak tidal strength of encounters is plotted against 
mass of the harassing galaxy. There is no indication of stronger 
tidal encounters occurring at the high end of the galaxy mass 
function 



have this type of tidal history, so these can be considered 
as typical infalls. It should be noted that even tidal histo- 
ries that include peak tidal forces as large as 3 Roche limits 
(occurring in two-thirds of infalls) do not show significantly 
stronger effects of harassment (see Table [2]). Significant tidal 
encounters that cause strong transformation or stellar losses 
occur for tidal forces that exceed ~ 4 Roche limits in the full 
numerical simulations. These occur in less than 25% of the 
infalls and cannot be considered statistically typical. Hence 
mild harassment can be considered the norm for newly ac- 
creted dwarf galaxies, at least for the orbits considered so 
far (see Section r6.4l for the dependency of this conclusion on 
orbital parameters). 



6.2 What drives strong tidal encounters? 

We are also in a position to answer the question of whether 
strong tidal forces are dominated by encounters with the 
rarer, more massive harassing galaxies. In Figure [S] the 
strength of the encounter versus the mass of the harassing 
galaxy is plotted. There is a large degree of scatter in the plot 
and little trend for stronger tidal encounters with increased 
mass of the harasser. Strong tides are not predominantly 
caused by encounters with the rarer, more massive cluster 
galaxies. A far stronger trend is seen when we instead plot 
tidal strength against the minimum separation between a 
harassing galaxy and the infalling dwarf - see Figure (9] This 
clearly indicates that strong tides are in fact driven by close 
encount ers with harasser ga laxies. This is contrary to state- 
ments in I Moore et al.l (|l996l ) suggesting that the bulk of the 
evolution is caused by encounters with galaxies brighter than 
~ L*. In fact the strongest encounters seen may not be well 
modelled by the simulations in this paper, as the proxim- 
ity of the galaxy-galaxy encounter is such that a collision of 
the baryonic contents is likely. If the colliding galaxies both 
contain a dissipative gas component, then there is the possi- 
bility for a dissipative (wet) merger. Fortunately, this is not 
a common event. < 6% of all encounters occur at separations 
of less than 15 kpc. The dwarf galaxy model initially has a 
gas disc that extends to ~ 15 kpc, hence passes closer than 
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Figure 9. Tidal strengtii of encounters is plotted against min- 
imum separation between tiie dwarf galaxy and the harassing 
galaxy. There is a strong trend for increasing strength of tidal 
encounters with decreasing separation. 



this, between similar galaxies, are likely to cause a collision 
rather than a tidal encounter. The maximum tidal strength 
encountered at a separation greater than 15 kpc is 5 times 
the Roche limit of the dwarf galaxy. Thus infalling galaxies 
who experience tidal forces as strong as this (e.g. Run 5) 
can still be considered to be well modelled by the simula- 
tions in this paper. However, it is not unreasonable for us 
to assume that the infalling dwarf galaxies are ram pressure 
stripped of their gas content before the collisions occur. In 
this case, the colhsionless stellar component of the dwarf 
galaxy will only be affected by gravitational tides, as mod- 
elled in these simulations, and the results of even stronger 
tides would still remain reasonable. This has interesting im- 
plications for the possibility of colhsionless (dry) mergers. 
Newly infalling galaxies are not captured by the tidal po- 
tentials of the harassers galaxies, but merely deviate from 
their original trajectory. 



6.3 Location and speed of encounters 

While ram pressure stripping is believed to be only signif- 
icant in the inner regions of the cluster, how common are 
high speed tidal encounters in the outskirts of the cluster? 
In this model, harassment is found to be similarly preferen- 
tial to the cluster centres. Statistically, the average radius 
of an encounter is ~ 200 ± 150 kpc within the cluster. Only 
~ 7% occur at radii beyond 500 k pc (beyond the core radius 
Eis defined in ISabatini etaP (|2005l );i. Hence, the effects of ha- 
rassment are very unlikely to be observed in the outskirts of 
the cluster, unless the dwarf galaxy has already made a pass 
of the cluster ce ntre (i.e. it is a member of the back-splash 
l|Gill et al.ll2005h galaxy population). 

The average velocity of an encounter is very high with 
a large standard deviation: ~ 2000 ± 600 km s~^. Hence 
the short-lived strong tides seen in Figure [2] are virtually 
the standard and consequently galaxy-galaxy capture events 
are very rare . Chance low speed encounters are vanishingly 
small, with < 4% occurring at velocities of < 500 km s~^. 
At first glance, these encounters appear to occur at very 



high velocity, considering the total velocity dispersion of the 
cluster is ~ 1000 km s~^. However, the elliptical orbits of 
the harasser galaxies causes the velocity dispersion in the 
central region of the cluster to be higher and this is where 
the majority of encounters occur. In addition, the plunging 
orbit of the infalling dwarf is far from the orbit of a virialised 
galaxy population raising encounter velocities further. 



6.4 Strong sensitivity to orbit 

6.4-1 A less plunging orbit 

The orbit utilised in the full numerical simulations and in the 
Monte-Carlo simulation described above, is not un-realistic 
but does involve a plunging orbit. We repeat the Monte- 
Carlo simulation but this time with a less plunging orbit. 
This allows us to test how strongly the tidal histories of 
infalling dwarf galaxies depend on their orbit through the 
cluster. 

For this we use the average orbital parameters for sub- 
str ucture in cluster -mass objects of a ACDM simulation - 
see et al. I l|2004l l. Despite a wide variety of cluster prop- 
erties, the orbits of the satellite populations were found to 
be similar with an average ellipticity e = 0.61, and an av- 
erage pericenter distance Vperi = 0.35r„ir- (with minimal 
scatter in both quantities). Here ellipticity is defined as 



e = 1 



peri / ' apo 



and 



is the virial radius of the clus- 



ter. For the cluster model presented in this paper, this fixes 
rperi = 385 kpc, and Vapo = 987 kpc. Hence the orbit is 
less eccentric, and does not plunge so deeply into the cluster 
centre. To produce this orbit, all of the Monte-Carlo parti- 
cles representing galaxies are positioned at the apocenter of 
the orbit, and given a 450 km s^^ azimuthal velocity. This 
produces a peak velocity at pericluster distance of ~ 1250 
km s~^. 

The result is significantly less strong harassment. In Fig- 
ure [To] a cummulative frequency plot of peak tidal forces 
encountered by a dwarf galaxy following a typical orbit is 
shown. Now 70% of galaxy infalls never experience a tidal 
force greater than two times the Roche limit, compared to 
50% for the previous more plunging orbit. Such a tidal his- 
tory typically resulted in only minor stellar losses, and minor 
disc heating in the full numerical simulations - see Table 
[2] Still Run 4 and Run 7 represent the most statistically 
common results of harassment. Strong tidal encounters that 
cause significant stripping and morphological transforma- 
tion (at tidal forces ~ 4 times the dwarf galaxies Roche limit 
or more) occur in less than 15% of the infalls, compared to 
25% for the previous plunging orbits utilised. Clearly, the 
infiuences of harassment are strongly dependent on the or- 
bital parameters of the orbiting galaxies, and this can be 
simply understood - the deeper the orbit plunges in to the 
cluster, the more harassing galaxies are encountered. How- 
ever, there must be an additional dependency on the time 
spent in the cluster core also. Hence it might be expected 
that the strongest harassment will occur for much less ec- 
centric orbits, that are close to circular, and remain near to 
the cluster center throughout their orbits. 

This does indeed seem to be the case for th e 
dwarf-harassment simulations in lMastropietro et al.l l|2005h . 
Galaxies that have large apocentric distances (~ r„ir) all 
only suffer very mild harassment with similar dark matter 
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Figure 10. Cummulative frequency plot of peak strength tidal 
encounters experienced by infalling dwarf galaxies on typical or- 
bits for substructure within ACDM simulations of cluster mass 
dark matter haloes 
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Figure 11. Cummulative frequency plot of peak strength tidal 
encounters experienced by infalling dwarf galaxies initially on a 
circular orbit, close to the cluster centre at radius 200 kpc 



and stellar losses to the dwarf galaxies in these simulations 
(~ 10% stellar and ~ 60 - 70% dark matter lost). Although 
the harassment model used in this study involves analytical 
calculation of tidal forces, it produces similar results to a full 
ACDM simulation of a cluster, although a direct compari- 
son is not possible due to the significant ly lower mass of the 
dwarf models in this study. Meanwhile in Mastropietro et al] 
almost all orbits whose apocentre is ~ 0.1 r„ir, and 
hence never leave the dense cluster core suffer significant 
mass loss (~ 50% and greater stellar loss, and ~ 95% dark 
matter loss). A low eccentricity, small apo-centric orbit is 
additionally more likely to make a number of repeat passes 
of the cluster core. If a galaxy makes a second pass, the 
chances of a strong encounter should double if the clus- 
ter mass and number of harassers has not changed signif- 
ic antly over this time. The s trongly harassed dwarf galaxies 
in lMastropietro et al.l |200^ complete numerous orbits over 
the lengt h of the simulati on. Similar conclusions are also 
drawn in iGill et al.l (|2004l ). However, this form of orbit is 
highly unlikely to occur for galaxies that are newly accreted 
into the cluster at the current epoch. The cluster now has 
a larger mass, and the galaxy's orbit must have high eccen- 
tricity to pass close to the cluster core where the effects of 
harassment are strong. This results in a long period orbit 
with a large apo-centre, where it is unlikely that the galaxy 
will make a second pass in > 6 Gyrs (for the NFW back- 
ground potential used in this study). Dwarfs that make nu- 
merous repeat cluster center passes must have formed close 
to the cluster core, or at least infallen when the cluster was 
significantly less massive. Such dwarfs could form their own 
sub-class of dwarf, as discussed in the 'Discussion and Sum- 
mary' section. 



6.4.2 A small apo-centric orbit 

To test the influence of harassment on galaxies that are not 
newly infalling into the cluster, we repeat the Monte-Carlo 
simulations only this time we place each Monte-Carlo parti- 
cle closer to the centre of the cluster, at a radius of 200 kpc. 



Run 


Apo-Peri 


DM 


Star 


{Vpeak/cr) 


(M/L)dy„ 


Run 1 


220-100 


12 


64 


38 


19 


Run 2 


200-50 


1 


3 






Run 3 


240-70 


1 


8 






Run 4 


240-120 


28 


85 


55 


33 



Table 4. Summary of effects of a small apocentric orbit (r<250 
kpc) harassment on the standard dwarf galaxy model; column 1 
is run description, column 2 is apocentre and pericentre of the 
orbit, column 3 and column 4 are percentage of dark matter and 
stellar mass at the end of the simulation, coiiimrt 5 and column 
6 is final (V^eafc/o")! and final {M/L)d;/n as a percentage of the 
orignal 



Each is provided with a purely azimuthal velocity of 660 
km s^^ producing an orbit that, in the absence of harasser 
galaxies, is virtually circular (zero eccentricity). A cummu- 
lative frequency plot of the peak tidal force experienced by 
each of the Monte-Carlo particles is presented in Figure [TT] 

Clearly, strong tidal encounters occur significantly more 
often for such a small apo-centric orbit. There are now no 
Monte-carlo particles that experience a peak tidal force of 
less than one Roche limit over 2.5 Gyrs of evolution. In fact 
>70% of all particles experience at least one tidal encounter 
of force greater than 4 times the Roche limit of the dwarf, 
which has been shown to cause significant dark matter and 
stellar losses. Hence for this type of orbit, strong harassment 
is the norm. To test the effects of such tidal encounters on 
small dwarf galaxies, we place the standard model dwarf 
galaxy at a 200 kpc radius and provide it with the same 
660 km s~^ aximuthal velocity. We then evolve the dwarf 
galaxy for 2.5 Gyrs under the infiuences of harassment. We 
repeat this test 4 times, with a different velocity vector each 
time (although all have a purely azimuthal velocity of 660 
km s~^) so as the 4 models each under go different tidal 
histories. The results can be seen in Table 3] 

Run 1 suffers strong harassment, in a similar manner 
to the rare strong encounters of the more plunging orbits, 
(see Run 6 of Table [2]) with ~ 90% dark matter losses, and 
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a substantial decrease in (V^eafe/o") and (M/L)dan- Run 2 
and Run 3 are both destroyed by tidal encounters, making a 
measurement of their final iVpeak/cr) and (M/L)dyn impos- 
sible. The Monte-carlo tests demonstrate that disruption by 
the cluster background potential alone is a rare ocurrence. 
<1% of Monte-carlo particles experience a tidal force above 
90% of the models Roche limit from the cluster background 
potential alone. Run 2 and 3 are dismantled by repeated 
strong tidal encounters and pass the cluster centre several 
times over the period of their orbit. Run 4 suffers mild ha- 
rassment - similar to the typical harassment of the plunging 
orbits (e.g. Run 1 and Run 7 of Table [2|. It can be consid- 
ered as one of the <30% of the Monte-carlo particles that 
avoids tidal encounters of greater than four times the Roche 
limit, hence is fairly atypical of harassment in this regime. 

In general, the dwarfs are typically destroyed or suffer 
significant dark matter and stellar losses when undergoing 
small apocentric orbits. The final stellar disks of the two 
surviving galaxies show significant heating, and appear to 
have increased in vertical scale-height substantially. It is es- 
timated that, if viewed edge-on, the ratio of their short-axis 
to long-axis would be ~ 0.8 and ~ 0.5 in Run 1 and Run 4 
respectively - hence Run 1 has become close to spherical in 
it's stellar distribution. 

Once more, a direct comparison between the results of 
iMast ropietro et al.l (|2005l ) is not possible due to the differ- 
ing mass-regime of the dwarf galaxy models, but substan- 
tial dark matter and stellar losses of equivalent magnitudes 
(~90 % and ~50 % respectively, and gr eater) are seen. For 
compa rison Gal 1,2,4,7-12 of Table 1 in IMastropietro et al.l 
(2005) have approximately similar orbital characteristics 
and mass-loss. However, one difference is that none of Mas- 
tropietro's dwarfs were considered completely destroyed by 
the end of their simulation. The lower mass dwarfs are more 
susceptible to disruption by the combined action of the 
background cluster potential (which is continuously stronger 
than in the plunging orbit case), and numerous strong tidal 
encounters. In this situation, dark matter that is liberated 
by a strong encounter, is quickly swept away by the stronger 
tidal field of the background cluster potential. This brings 
into question the fate of the dark matter and stars from 
such disrupted galaxies. One possibility is a fraction of the 
dark matter and stars is absorbed into the central ellipti- 
cal of the clusters, while any stars that are not absorbed in 
this manner contribute to the l ow surface brightn ess back- 
ground observed in the cluster (|Mihos et al.ll2005l ). Indeed, 
the filaments and streams of stars observed at low surface 
brightnesses between the giant central elliptical galaxies in 
Virgo, may have their origin in the merging of heavily dis- 
rupted harassed dwarfs with the giant ellipticals. 



7 DISCUSSION AND SUMMARY 



As discussed in lLisker et al.l (|2007l ). it is important that we 
recognise the various sub-classes of dwarf galaxy if we are 
to understand their origin. If we split the cluster dwarfs into 
nucleated and non-nucleated dwarfs, we find significant vari- 
ations between the two populations. In order to recover the 
intrinsic, three dimensional shape of dwarf galaxies, fiatten- 
ing distribut ions combined with the su r face brightness tes t 
are utilised l|BinggeH fc Popescul [l995l . iLisker et al.ll2007l) . 



Comparisons of non-nucleated dwarf ellipticals, with dwarf 
irregulars and late-type spirals show no statistically sig- 
nificant difference in their flattening distributions as mea- 
sured by a standard Ko lmogorov-Smirnov test. However 
iBinggeh fc Popescul (|l995l ) find that nucleated dwarfs m aybe 
as oblate as classical ellipticals, and lLisker et al.l l|2007l ) find 
they are more akin to E3 or E4 type giant ellipticals. Fur- 
ther evidence for a different evolutionary history is reflected 
in their spatial distribution. Nucleated dwarfs are found to 
be concentrated towards the centre of the Virgo and Fornax 
cluster like the giant E and SO galaxies. Meanwhile non- 
nucleated dwarfs are distributed like the late-type spirals 
and irregulars, sug gesting they are far fr om virialised within 
the cluster (iFerguson fc Sandagd 19891 ). A similar conclu- 
sion is drawn in iLisker et all (|2007l ~ who uses radial veloc- 
ity proflles to suggest that the nucleated dwarf ellipticals 
are indeed virialised, whereas the non-nucleated dwarfs still 
bear the signatures of a recent infall. 

This accumulation of evidence could suggest that nu- 
cleated dwarfs formed in situ within the cluster, whereas 
the non-nucleated population have been accreted from the 
cluster outskirts. The research presented in this paper sug- 
gests that harassment does not provide a significant mecha- 
nism to convert newly infalling late-type dwarfs into dwarf 
ellipticals. Ram pressure stripping of late-type dwarfs may 
provide a far more effective mechanism to form the non- 
nucleated dwarf s (.Boselli et al. 2008 ) and will be studied in 
further detail in lSmith et al.l I 201ol ) in combination with the 
harassment model presented in this paper. 

However, it is less clear how ram pressure stripping 
could have formed the nuclei fo und in the nucleated dwar f 
population. The simulations of IMastropietro et al.l (|2005l ) 
demonstrate that harassment may be far more effective on 
a more centrally concentrated, and virialised dwarf popu- 
lation that formed in situ in the cluster. The early cluster 
environment would be expected to have lower intra-cluster 
medium densities, and so ram pressure stripping would be 
far less efficient than it is today. If their inter-stellar medium 
is not stripped, then ti dal encounters c ould efficiently drive 
gas into their centres (Mo ore et al. _1998i), and could facili- 
tate the formation of central nuclei. Perhaps these galaxies 
were the progenitors of current day nucleated dwarf ellipti- 
cals whose oblate shapes may refiect the stronger harassment 
they have incurred. 



7.1 Key results 

(i) In general (> 75% of infalls), dwarf galaxy models 
that infall from the outskirts of the cluster suffer only mild 
harassment. Typical stellar losses are ~ 10%, while typi- 
cal dark matter losses are ~ 60%, resulting in a reduction 
in dynamical mass-to-light ratio by a factor ~ 2. Typical 
(Vpeafc/o") ratios are also reduced by approximately the same 
factor, but final stellar discs remain rotationally dominated, 
and appear as thickened discs, and often with tidally induced 
spiral structure. 

(ii) In rarer cases (< 25% of infalls), dwarf galaxy models 
infalling from the cluster outskirts suffer strong harassment, 
losing ~ 70% of their stars and ~ 90% of their dark matter. 

(iii) A Monte-Carlo simulation reveals that strong harass- 
ment occurs in < 25% of infalls for the deeply plunging orbit, 
and < 15% of infalls for a more typical orbit of substructure 
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within a ACDM cluster-mass halo. The effects of harass- 
ment are strongly dependent on orbital parameters, as this 
dictates time spent in the cluster centre and the number- 
density of harassers encountered. 

(iv) For a low eccentricity, small apo-centric orbit, strong 
harassment is typical, occurring in > 70% of galaxies on 
such orbits over 2.5 Gyrs. This typically results in strong 
mass loss (> 90% dark matter, > 50% stars) or complete 
disruption of these low mass dwarf models. 

(v) Loss of stellar mass is caused by impulsive heating of 
the stellar disc by short-lived tidal interactions. Stars that 
are stripped from the disc by harassment form very low sur- 
face brightness streams {^b > 31 mag arcsec"'^) 

(vi) The Monte-Carlo simulation reveals that average 
tidal encounters occur at high velocities (~ 2000 km s~^), 
with chance low speed encounters being very rare (< 4% of 
encounters have a relative velocity < 500 km s~^). 

(vii) For a plunging infall, galaxy-galaxy encounters are 
centrally concentrated, typically occurring at a radius ~ 
l/5th the cluster virial radius. Encounters in the outer clus- 
ter are rare (~ 7% occur at a cluster radius > 500 kpc). 
Therefore the effects of harassment are unlikely to be ob- 
served in the outskirts of the cluster for newly accreted 
galaxies - with the exception of the back-splash population 
that has already completed one near pass of the cluster cen- 
tre. 

(viii) Strong tidal interactions are not primarily caused 
by encounters with massive harasser galaxies, but are driven 
by near-pass encounters with more typical cluster galaxies. 

(ix) Spira l structure such as those observed in 
iLisker et all (2006) can be induced by tidal encoun- 
ters, and therefore is not necessarily an indication that that 
original galaxy contained the same spiral structure. 



8 CONCLUSIONS 

Harassment of small late-type dwarf galaxies that have re- 
cently been accreted mto the cluster is not an effective mecha- 
nism for forming the thickened discs of dwarf ellipticals. De- 
spite their low mass, they are additionally very concentrated 
in mass, and hence are ro bust against the ef fects of harass- 
ment. In previous studies ()Moore et al.lll99^ . harassment is 
the effect of repeated long-range fast encounters. For dwarfs 
in the low-mass but more concentrated regime studied here, 
harassment is only signicant in repeated, or even individual, 
close-range high speed encounters where tidal forces are the 
strongest, and these are rare for objects infalling into the 
cluster at the current epoch. Currently it would be a chal- 
lenge to observe the typical effects of harassment on dwarfs 
that have recently infallen into the cluster (very low surface 
brightness extragalactic features, mild disc thickening, and 
a ~ 50% reduction in dynamical mass-to-light ratios) . How- 
ever harassment appears far more effective on dwarf galax- 
ies that have formed in situ within the cluster as it formed, 
possibly forming nucleated cluster dwarfs. Alternatively ha- 
rassment of larger, low surface brightness spirals into dwarf 
ellipticals in the manner suggested bv Moore et al. (1998) is 
another possibility. However, the simulations of lMoore et al.l 
l|l999l ) demonstrate the harassment is not effective on all 
disc galaxies - while a giant low surface brightness disc is 
effectively tranformed, a high surface brightness giant like 



our own galaxy is only mildly affected. Therefore if cluster 
dwarf ellipticals have their origins in larger disc galaxies, 
a very large population of specifically low surface bright- 
ness field disc galaxies is required, in order to form the vast 
numbers of dwarf ellipticals seen in clusters today. The for- 
mation of dwarf ellipticals from dwarf irregulars via ram 
pressure stripping is another convincing possibility as dis- 
cussed in lBoseUi et~aLl ([2008'). Ram-pressure stripping need 
only terminate their star formation to produ ce objects that 
appe ar very much like dwarf ellipticals today l|van Zee et al.l 
|2004 ). 
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